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Abstract 

We analyze the decay modes B° D*-D°K + and B° D*-D*°K + and, 
using the available experimental data, we find bounds for the constants g and 
h describing the strong coupling of heavy mesons to light pseudoscalar mesons. 
Both the decay channels are dominated by broad L = 1 charm resonances; 
the dominance is effective also in B° — > D~D°K + and B° — > D~D*°K + . 



1 Introduction 



Recently, the BaBar Collaboration has observed the three-body B decay modes 

D*-D°K + (1) 
D *- D *o K + ( 2 ) 

= (2.8 ±0.7 ±0.5) HT 3 (3) 
= (6.8 ± 1.7 ± 1.7) HT 3 . (4) 

a measurement has also been reported by the 

= (3.2 ±0.8 ±0.7) 10~ 3 . (5) 

The results ([!]), (|]) and (0) (although preliminary) represent a significant improvement 
with respect to the previously available data, obtained by the CLEO Collaboration: 
B(B° -> D*-D°K+) = (0.45lg^ ± 0.08) 10~ 2 and B(B° -> D*-D*°K + ) = (1.30jl[j$ ± 
0.27) 10~ 2 ||. We expect that the experimental analysis of the processes (|1|) and (0), 
including the study of the Dalitz plot, will be further pursued in the forthcoming future. 
The purpose of this note is to interpret the observations made so far, using the results 
(H)"@ together with the existing datum on the two-body decay B° — > D*~D* + Q: 

B(B° -> D*~D* + ) = (19 ±6) 10" 3 . (6) 

The interest for B transitions into a pair of D^*> and a Kaon is manifold. It has been 
proposed to use the modes B°(B°) -> D-D + K s and B°(B°) -> D*-D* + K s (analogous 
to (|D) to investigate CP violation effects in neutral B decays at the B factories ||. Such 
processes are induced at the quark level by the transitions b — > ccs and b ccs and are 
Cabibbo-favoured as in the case of B°(B°) — > J/ipKg, with a tiny penguin contribution. 
Studies of the time-dependent Dalitz plot would provide us with information about the 
weak mixing angles, namely the phase (3 related to the B° — B° mixing. In particular, 
since amplitudes with different strong phases corresponding to various intermediate states 
contribute to the three-body B — > D^'D^Kg decays, one envisages the possibility of 
measuring both sin(2/3) and cos(2/3) by suitable Dalitz plot analyses 0. 



B° - 
B° - 

measuring the branching fractions Q: 

B{B° -> D*-D°K + ) 
B(B° -> D*-D*°K + ) 

For the decay channel B° —> D*~D°K^ 
Belle Collaboration 0: 

B{B° -> D*-D°K + ) 
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Another reason of interest concerns the possibility of carrying out tests of factorization 
for nonleptonic B decays. It is reasonable to assume that the modes (|IJ) and (Q) mainly 
proceed through two-body intermediate states, such as 

B°^D*-D? , (7) 

followed by the strong transition 

D x _ D (*)o K +. ( 8 ) 

Df are charmed strange mesons; a typical diagram is depicted in fig.|l[ In the factorization 
approximation the amplitude of the process in (0) is expressed as the product of the 
semileptonic B° — > D*~ matrix element and the Df current- vacuum matrix element. 
In the infinite charm quark mass limit, the only contributions with non- vanishing D* 



£)*- £)(*)0 




Figure 1: Diagram contributing to the decay B° — > D* D ( -*'°K + . The box represents a 
weak transition, the dot a strong vertex. 



current-vacuum matrix elements correspond to the states = D* and D sQ (with their 
radial excitations) for B° — > D*~D°K + , and Df = D*, D s and D* sl (together with 
their radial excitations) for B° — > D*~ D*° K + . D s q and D* sl are positive parity mesons 
belonging to the sf = | + heavy meson (sc) doublet, sf being the spin-parity of the 
light degrees of freedom in the meson. Therefore, the number of independent amplitudes 
contributing to ([]]) and (fj) is limited, and it is possible to study relations, e.g., with the 
mode B° — > D*~D* + for which the experimental datum (|6|) is available. Moreover, one 
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can exploit a possible dominance of the positive parity intermediate states to study the 
features of these (so far) unobserved states. 

There is a further reason of interest in the processes (|l|) and (|2]). If the main mech- 
anism for the three-body D*~ K + final states is the production of a pair of D*~ Df 
mesons followed by the strong transition (Rl), one can use such decay modes to access 
the couplings of heavy mesons to light pseudoscalar states. For these quantities little 
experimental information is currently available. The CLEO Collaboration has provided 
the first determination of the strong coupling constant go*Dn governing the transition 
D* + — > D°7T + , using the recent measurement of the total width of the D* + meson || 

r(ZT+) = 96±4±22 KeV 1 (9) 

together with the experimental branching fraction B(D* + — > D°tt + ) = (67.7 ± 0.5) 10~ 2 
Q. The result for the coupling, defined by the matrix element 

< D°(k)n + (q)\D* + (p, e) >= g D * Dlx e ■ q (10) 
(e is the D* polarization vector), is: 

9d*d-it = 17.9 ± 0.3 ± 1.9 . (11) 
Rewriting Qd*d-k m terms of an effective coupling g D : 

9D*D* = ^— f 9o = ^— f g(l + 0{ — )) , (12) 

one translates the result QTTj ) into 

g D = 0.59 ±0.01 ±0.07. (13) 

In the heavy quark limit the parameter g in (|1^) describes the strong coupling of charmed 
mesons as well as of beauty mesons to the members of the octet of light pseudoscalars; 
therefore, neglecting SU(3)f breaking effects, this parameter enters some matrix elements 
governing the transitions in (|8]). In addition, together with analogous couplings, g repre- 
sents a basic quantity in the heavy-quark chiral effective theory || [|, and therefore it is 
worth searching information about it from all available experimental data, and comparing 
the results with the predictions that vary in the rather wide range 0.2 < g < 0.7 [ID|. 
This is a purpose of the present note. 

In the next section we analyze the decay modes ([!]) and @ and discuss how to access 
the relevant strong couplings. Numerical results follow in section 3. The conclusions are 
drawn at the end of the note. 



1 This result updates the upper bound provided by the ACCMOR Collaboration: T(D* + ) < 131 KeV 
at 90% c.l.|j7). 



2 Decay modes B° D*~D°K + and £° D*~D*°K + 

Let us consider the processes ([!]) and (0): 

S°(p) - D*-(p_,e_)L> ( Pl? )^ + (g) 
S o (p) ^ D*~(p-,e-)D*°(p D *,e)K + (q) 

with the momenta p = m^f, p_ = m B *V- and £>£>(*) = m D i*)W expressed in terms of 
the heavy-meson four-velocities v, t>_ and w. Neglecting penguin contributions, the 
processes are governed by the effective weak Hamiltonian: 

H w = |u;«iM! -l*)ccf(l -75)5 (14) 

where GV is the Fermi constant, Vij are CKM matrix elements and the parameter a± reads 

a\ = (ci + with 01,2 short-distance Wilson coefficients and iV c the number of colors. 

\ i v c / 

Dalitz plot variables of the decays can be defined: 

s = (p D w + 

= (p_ + g) 2 (15) 

and a set of invariant variables, in terms of the four-velocities v,V- and w, can be intro- 
duced: 

m 2 B + m 2 D * — s 



v ■ v. 



v ■ w 



2m B m B) * 

m% + m 2 DM - s 

2m B m DM 

m 2 B + m 2 K - s - s_ 
f _ ■ w = (16) 

2m D *m D (,) 



v ■ q 
V- ■ q 
w ■ q 



s + s_ - m D * - mj, w 
s_ — m|)» — m 2 ^ 
s - m\ - m D (») 



2m D M 

In the plane (s, s_) the accessible kinematical region is defined by the conditions 

(m D (*) + m K ) 2 < s < (m B - m D *) 2 

(0-< s <(«-)+ ( 17 ) 

5 



where 

(s_)± = m D * +m? K - — 



s-m B + m D ,){s + m K - m DMj 



2s 

^A 1/2 (s, m 2 K , m 2 DM )X 1/2 (s } m 2 B , m 2 D *) 

A being the triangular function. 

We assume that the decays B° — > D*~ K + proceed through polar diagrams such 
as the one depicted in fig.|l], computed adopting the factorization approximation for the 
weak transition. In the case of B° — > D*~ D°K + , the pole can be either a vector ( J p = 1~) 
meson: D* + , or a scalar (J p = + ) meson: Df , with their radial excitations. For the 
decay B° -> D*-D*°K + , the possible poles are: D* + (J p = 1"), L>+ (J p = 0~) and D\t 
(J p = 1 + ) and their radial excitations. Therefore, the calculation of the amplitudes in 
fig{l] requires the strong vertices 

<D\p D )K + (q)\D:(p Dt ,e s )> = g D * sDK (e s ■ q) 
<D°{p D )K + (q)\D+(p Ds0 )> = g DsoDK 
<D*°(p D *,e)K+(q)\D;(p Dt ,e s )> = i 9 -^f- e^e^pt^ (19) 

<D*°(p* D ,e)K + (q)\D s (p Ds )> = g D . D . K (e* ■ q) 
< D*\p* D ,e)K + (q)\D:i(p D:i ,e s ) > = (e* • e s ){p D * ■ q) , 

si 

and analogous matrix elements involving radial Df resonances. In the heavy quark limit, 
all the couplings in (|19D can be expressed in terms of two different coupling constants 
g and h, for negative and positive parity states, respectively. This can be shown 
considering the effective lagrangian describing the interactions of heavy mesons with the 
light pseudoscalars. In the limit mq — > oo the heavy quark in the heavy mesons only 
acts as a static colour source, and the gluons decouple from the heavy quark spin sq, 
thus implying a SU(2Nf) spin-flavour symmetry 0, [12|. At the opposite energy scale, 
for vanishing masses of the up, down and strange quarks, the QCD SU(3)l x SU(3)r 
chiral symmetry is spontaneously broken, the Goldstone bosons being the octet of the 
light pseudoscalar mesons. Both the heavy quark spin-flavour and the chiral symmetries 
can be realized in a QCD effective lagrangian ||, where the term describing the strong 
interactions of the heavy negative and positive parity mesons with the light pseudoscalars 
reads: 

d = % g Tr{H bW ALH a } + [ih Tr{H hlillh A^ a S a } + h.c. ] . (20) 
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The fields H a in fl20|) describe the negative parity J p = (0 ,1 ) qQ meson doublet, with 



[PLY - P«7fe] 



(21) 



2 L~ a/i 

the operators and P a respectively annihilating the 1~ and 0~ mesons of four-velocity 
v (a = u, d, s is a light flavour index). Analogously, the fields S a describe the positive 
parity states, with sj 



i + . 



(22) 



The octet of the light pseudoscalar mesons is included in (ED) through the field £ = e 
with 



V 



7T 

if" 



-A/ivr 



K c 



(23) 



I 7 ? / 



and f n = 131 Mel 7 . Finally, the operator A in (120) reads 



ba 



(24) 



From the definitions in ( B5 ) and from eq. ([20J it is straightforward to derive the relations 

2 y /m D .m D 



9d*dk 
9d s0 dk 
9d*d*k 
9d*d s k 
9d* s1 d*k 



K 



-^/m Ds0 m D 



m D s0 ~ m D h 



™d s0 f. 



K 



2m 



Dt 



f 9 
JK 

2^/m D *m Ds 



(25) 



A 



2^/m D *m D * 
Tk ' 

where we have kept some SU(3) flavor breaking terms in the masses of the mesons 
and in the leptonic constant Jk- 

On the other hand, in the factorization approximation the calculation of the weak 
transition ([?]) requires the semileptonic B° — > D*~ matrix element and the decays constant 
of the poles Df. In the heavy quark limit, the former is given in terms of the Isgur-Wise 
function £: 

<D*~{v^eJ)\V IJ --A li \B{v)> = ^/m B m D « £(v ■ vJ) ti a 

( -i6 pa ^v p v^x-{l + v-v^)g a ^ + v a vi) (26) 



while the decay constants are defined by 

< Dj( PDs )\c^(l- l5 )s\0> = if D y Ds 

<D* s + (p D: ,e s )\^(l-^)s\0> = f D: m D .e? 

<Dt ( PDs0 )\cY(l-i 5 )s\0> = if Ds y Ds0 (27) 

< J D: i + (^ i , es )|c 7 ^(l-75)s|0> = f Dt m Dti eT. 

In the heavy quark limit, the leptonic constants fo s and fo*, as well as fo s0 and Jd* s1 , are 
simply related. 

It is now straightforward to work out the amplitude of B° — > D*~D°K + proceeding 
via the D* intermediate state: 

iJC f D *m D * g D * DK t , N , oc x 

( _ M (m g t; - m D *v_) ■ q (m B v - m D *t;_)^ \ 
V 9 m Dt rn D * I 



{ie IIua f3V*v f3 - g^il + v • V-) + v v {v-)^ 



G F 

with /C = —=V C sV* h aiy/rriBrriD* and T^* the Z)* decay width. Analogously, the amplitude 
V2 

.4.2 relative to the D s0 contribution to B° — > D*~D°K + reads: 



^ fD s0 m Ds0 g Ds 



A 2 = ~J x»*dk £( v . v _ )e *_v ( 29) 

s - m D s0 + im Ds0 r Ds0 

f (m B u-m ]) .ti_L i , (m B t> - m D *v_) ■ v_ } 

< (1 + V-V-) + v u } . 

{ m D s{) m Ds0 J 

As for B° — > D*~D*°K + , the amplitudes ^4*, ^ and .4£ corresponding to the D*, D s and 
D*! intermediate states are given by: 

A * fD*m D * g D * D * K ^ „ , fl m D . . ^, 

s — TTijj* + mi£)*L vr mo* 

{i e II uai3V^v /3 - 5^(1 + v ■ v-) + ; (30) 

^2 = D ' K r £ (u • O ("tB + (e ■ g) (e* ■ ^) ; (31) 

s-m D +im Da T Da 



)CfD* 1 'mD*gD* 1 D*K 

s-m 2 Dti +im Dt r Dti 

{ - ie TUaP v"v p + g UT (l +vv-)- v v (v-) T 

(m B v - m D *V-) v (m B v - m D *V-) T 

(1 + V • V-) 

m D * m D * 

si si 



(m B v - m D *vJ) ■ v_ (m B v - m D *v^) T 



si si 



Expressions analogous to eqs.(p8|)-( p^) , with appropriate masses, widths, leptonic con- 
stants and strong couplings, hold for the contributions of the radial excitations of negative 
and positive parity mesons. Such contributions are suppressed by the small numerical 
values of the leptonic constants and of the effective couplings. This can be shown, for 
example, using the relativistic constituent quark model in ref.|13|, where one obtains: 
/d'//ds — 0-73, D' s being the first radial excitation of D s . In the same model, using 



the method described in |H| , one obtains: Qr>' a D*K I 9d s d*k — 0.32. Analogous reductions 
occur for positive parity states. A further suppression is due to the large decay width of 
the excited states. Therefore, one can conclude that the first radial excitations contribute 
to the amplitudes of the processes (|l|) and @ by less than 15% with respect to the con- 
tribution of the corresponding low-lying states, an uncertainty that can be included in 
the error affecting the effective couplings we are studying in this note. 

3 Numerical analysis and discussion 

On the basis of the above considerations, we write down the widths of the decay modes 
B° -> D*-D°K + and B° -> D*-D*°K + as follows: 

T(B° - D-DWK+) = f imB ~ mDt)2 ds f is - )+ ds _ (33) 

J{m DM +m K ) 2 J(s-)- dsdS- 



rp -i 

with -^^(B° -> D*~DM°K + ) = , ^ \ \AP and 
dsds_ K ' (27r) 3 32m|' 1 

A(B° -> D*-D°K + ) = A * 



i=l,2 

A(B° -> D*-D*°K + ) = ]T A* • (34) 

i=l,2,3 

The decay widths depend on the effective couplings g and h. They also depend on SM 
parameters, such as ^V cs V* b , on the leptonic constants fo a , . . ., on the Wilson coefficients 
c\ t 2 as well as on the Isgur-Wise form factor £. All such parameters appear in the same 
combination in the factorized amplitude of the two-body decay in (|6]): 

A ( B o D *- D * s ) = Kf D . mD: Z (v ■ v.) e* a eT[i e pct ^v x _v p + g afl (l +vv-)- v a v v \ . (35) 

Therefore, in the ratios 



B(B° -> /J*-/J* if +) 
^* " B(B° - D*~D* + ) (36) 

one gets rid of the dependences on V^, V^ s and ai. As for the Isgur-Wise function, the 
linear form 

£(vv-) = l-p 2 (v-v^ -1) (37) 

is well suited due to the narrow range of momentum transfer involved in the decays we 
are considering. A strong correlation has been observed between the measured values of 
Vd, and the slope parameter p 2 in the analyses of the semileptonic B° — > D*~£i> decay 
spectrum and in the studies of two-body B transitions in the factorization approximation 



151 |lq1 . However, in the ratios (|3q) such a correlation is essentially removed, and similar 

results are obtained varying p 2 in the range p 2 = 1.38 — 1.54. 

Sd Id Id* 

Concerning the leptonic constants, we use — s - = 1 and "° = sl = 1. The former 

Id* Jd* s Jd* 

ratio exactly holds in the infinite charm quark limit. The latter one allow us to reduce 

the number of input parameters, since a deviation from unity can be reabsorbed in the 

numerical result for the parameter h. For the masses of the excited charm mesons, we 

use rnn s0 = m D* sl = m D s ± A, with A = 0.5 GeV |T7J . 



The final set of input quantities involves the decay widths of the intermediate states. 
One can neglect the D s width (T(D 8 ) = 1.33 ± 0.03 10~ 9 MeV |§). As for D * s , as wel1 as 
for the positive parity charmed states, the widths depend on the effective couplings g and 
h. Using the experimental branching fractions B(D* + — ► Dfn ) and B(D* + — > Df^), 
together with the central value of g in (0), one obtains T(D*) = 1.03 MeV; consequently, 
we use the expression T(D*) = 1.03 (pfg) 2 MeV in the analysis for constraining the strong 
coupling g. Moreover, assuming that the decay widths of the positive parity states are 
saturated by two-body transitions, one gets T(D s q) = 180 (p^g) 2 MeV and T(D* 1 ) = 



165 (o^g) 2 MeV 



Solutions of the equations 



RD = R D \exp = 0.15 ±0.07 
R D *=R D *\exp = 0.36 ±0.17 (3* 



in the variables (g,h) are found, considering the central values in (|38|), for (g,h) = 
(0.05,-0.59) and (g,h) = (0.0, ±0.60). The solutions for g are smaller than the re- 
sult in (0), while the results for h are compatible with the theoretical expectations 
h = —0.52 ± 0.17 and h = —0.56 ± 0.28 |T7|]. However, before drawing conclusions from 
these results, it is worth analyzing the 1— and 2 — a regions in the plane (g, h), obtained 
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considering the experimental errors in flSSP - Such regions are depicted in fig.@. They are 
rather tightly bounded along the h direction, while the dependence on g is mild and the 
range of the allowed values of g extends over all the values between g = and the CLEO 
result eq.(p"3|). Along the h axis, the allowed regions correspond to \h\ = 0.6 ± 0.2. The 
conclusion is that the main contributions to the processes ([!]) and (f2|) are related not 
to the 0~ and 1~, D s and D* intermediate states, but to the positive parity + and 1 + 
states D s0 and D* l} since the amplitudes display a minor sensitivity to the coupling of the 
negative parity intermediate states. This is interesting from the phenomenological point 

l 

0.5 
43 

-0.5 
-1 

0.2 0.4 0.6 0.8 
g 

Figure 2: 1 — o (continuous lines) and 2 — o (dashed lines) region in the (g, h) plane, as 
obtained from the ratios Rd and Rd* in (|36|) . The vertical lines represent the result (|i~3|); 
the shaded area corresponds to the region excluded by the upper bound g < 0.76 from 
ref.0. 



of view, since it implies that three-body B° — > D*~D°K + and B° — > D*~D*°K + decay 
modes are well suited for separately studying the properties of the (so far unobserved) 
D s q and D* sl resonances ||. The analysis can be done by studying the Dalitz plot of the 
three-body decay. For the mode B° — > D*~D°K + the expected differential decay width 
is depicted in fig.^. It has been obtained for g = 0.5, h = —0.6, a\ = 1.1, together with 
V cb = 0.04 and V cs = 0.974 |§. As for f Ds , we use the value f Ds = 240 MeV obtained 
from the fit of (||); it is compatible, within the errors, with the value reported by ||: 
fo s = 280 ± 19 ± 28 ± 34 MeV. The distribution for g = 0.3 is completely similar. 

One notices that the main variation of the differential decay distribution occurs along 
the direction of the invariant D°K + mass, a feature related to the unique topology of 



11 



the (Cabibbo and color allowed) amplitudes governing the mode (|l|). The Dalitz plot for 
B° — > D*~ D*°K + , depicted in fig||, shows similar features. 




Figure 3: Differential decay width — - — (left) and Dalitz plot (right) of the decay 
B° -> D*-D°K+. Units of s and s_ are GeV 2 . 




s 



Figure 4: Differential decay width (left) and Dalitz plot (right) of the decay -B — > 
D*~D*°K + . Units are as in figure ||. 

The prominent role of the intermediate states D s q and D* sl makes the processes ([!]) 
and (0) particularly promising for the analysis of broad orbital excitations of the cs meson 
system. It is worth noticing, however, that also other three-body B decays can be well 
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suited for such an investigation. Examples are B° — > D~D°K + and B° — > D~ D*° K + , for 
which no experimental results are available. Since the matrix element of B° — > D~ can be 
related to B° — > D*~ in the heavy quark limit, it is possible to determine, in the scheme 
described in the previous section, the properties of the channels B° — > D~ K + . One 
predicts: 



B(B° -> D~D°K+ 
B(B° -»• D*-D°K~ 
8{B° -> D-D*°K+ 



2.11 

0.27 (39) 



-6(5° -> D*-D*°K+) 
which imply, considering the experimental data in (§)-(§): 

B(B°^D-D°K + ) = (6.3 ± 1.8) 10" 3 

B(B° -> D-D*°K + ) = (1.8 ±0.7) 1(T 3 , (40) 

in a range accessible to current experiments. The expected decay distributions and the 

Dalitz plots, depicted in fig.|5|, are similar to those of the modes ([[]) and (0), with features 

that it will be interesting to experimentally investigate. 

We conclude this section with a comment on the two main theoretical uncertainties in 

our approach, the use of the heavy quark limit both for beauty and charm quarks, and 

the factorization assumed for the nonleptonic matrix elements. The two uncertainties 

are correlated, and a quantitative assessment of their role is not a trivial task. If we 

consider them separately, we can presume that several corrections are compensated 

m Q 

in the ratios used as the basis of our analysis. As for factorization, the matrix elements 
governing the decays considered in this note are different from the matrix elements for 



which factorization has been proved in the infinite b mass limit ||18|| . Nevertheless, the 
study of various processes of the type considered here, i.e. color allowed B transitions 
to charm mesons, shows that factorization reproduces the available data within their 



current errors fllB] . Our analysis can be considered as a further test of factorization; 
experimental measurements will be helpful in shedding light on the size and the type of 
possible deviations. 



4 Conclusions 

The analysis of the decays B° — > D*~D°K + and B° — > D*~D*°K + shows that they mainly 
proceed through positive parity intermediate states, D s0 and D*±. Other contributions are 
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s 




D-D°K+ (up) and B° -> D-D*°K + (down). Units are as in fig.0. 



less significant, so that such three-body B° transitions appear to be well suited for studying 
the features of the low-lying orbital excitations of the (sc) meson system. Currently 
available experimental data allow us to constrain the strong coupling between such orbital 
excitations, the negative parity charmed mesons and the Kaon, in the region \h\ = 0.6±0.2, 
close to the expected values. An improvement in the accuracy of the measurements would 
further constrain this parameter. On the other hand, the coupling g is found in a range 
which extends from zero to the CLEO measurement ( |T3| ) . 

The Dalitz plots relative to B° -> D*-D°K + and B° -> D*-D*°K + are expected to 
display peculiar features, namely the main dependence on the invariant D^°K + mass, 
that can be experimentally tested. The investigations of other modes, B° — > D~D°K + 
and B° — > D~ D*°K + , can also provide us with information about positive parity charm 
states; their expected decay rates are accessible to current experiments. 
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